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V  -0E0  TRACKER  MODELS 

This  preliminary  report  considers  modem ng  of  two  classes  of  video 
trackers  from  their  algorithmic  and  measurement  aspects.  Because  they  are 
widely  used  and  easily  implementable  in  hardware,  the  classes  chosen  were  the 
binary  centroid  tracker  and  the  binary  correlation  tracker.  Many  variations 
exist  even  in  these  classes;  their  value  in  the  present  setting  is  to 
demonstrate  the  modelling  aspects  and  the  salient  features  that  may  be  derived 
from  an  analysis  of  the  behavior  of  these  systems  iri  the  presence  of  noise. 

The  models  predict  tracker  offset  errors  (mean  error)  and  track  point 
"jitter1*  (the  standard  deviation  of  the  error). 

Track  performance,  in  general,  is  dependent  on  the  character  of  the 
object  pattern  being  tracked  anc!  on  the  distribution  of  the  noise  associated 
in  the  formation  of  the  array  of  pixels  which  represent  the  pattern.  In  the 
analyses  which  follow,  a  simple  square  pattern  is  used  and  Gaussian  additive 
noise  is  assumed.  The  models  permit  any  geometric  patter;!  as  input  and  the 
question  of  noise  distribution  is  answered  by  the  specification  of  the 
probabilities  of  detection  and  false  alarm  in  the  formation  of  the  binary 
pattern  from  the  linear  data  input.  .  I  -  .. .  ,  t  .  ,-/v ,  %  r  , 

Data  is  presented  scaled  by  the  pixel  dimensions.  Additional  errors 
occur  due  to  spatial  quantization  which,  based  on  the  assumption  of 
uniformity,  adds  a  variance  of  1/1?  pixel  to  each  computed  answer.  This 
factor  is  included  in  the  curves  of  rms  error  versus  SNR;  further,  the  curve 
data  is  the  root  sum  square  of  the  individual  channel  errors. 

The  salient  features  of  the  analyses  are  as  follows: 

o  For  the  binary  centroid  tracker,  the  bias  error  is  proportional  to 
the  track  error. 

o  For  the  binary  correlation  tracker,  the  bias  error  is  a  function  of 
the  correlator  search  pattern. 

o  The  correlation  tracker  is  generally  better  than  the  centroid 
tracker.  A  curve  of  the  RSS  errors  is  shown  in  Figure  1. 


BINARY  CENTROID  TRACK F p 

The  analysis  that  follows  provides  an  estimate  of  the  error  in  the 
centroid  measuremen  of  a  target  pattern  due  to  the  presence  of  additive 
noise.  Certain  approximations  and  assumptions  are  evoked  in  order  to  ensure 
tractability  of  the  mathematics. 

The  input  to  the  tracker  is  assumed  to  be  provided  by  a  thresholded 
n  by  n  matrix  P  of  pixels  derived  from  the  imaging  sensor  output  data.  For  a 
threshold  level,  T,  the  binary  tracking  pattern  is  defined  as  follows: 

r  •  -  T  ■  i  ,J:-1»  •  •  .n 

1 J  -  -  m 


85 


The  "one"  values  of  the  pattern  resul t  from  two  sources: 

1)  Target  pattern  detections  with  probability,  Pf;  and, 

2)  False  alarm  or  background  detections  with  probability,  Pp. 

Row  and  column  vectors  are  formed  from  the  binary  matrix  to  compute  the 
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The  values  Xi,  y-f  t  and  S  are  random  variables  which  are  assumed  to  result 
from  a  series  of  independent  Bernoulli  trials.  The  number  of  "ones"  that 
occur  (or  the  integer  value  given)  is  described  by  a  binomial  probability  law: 

Probability  [n=k]  =  (  p*  (l-p)01-*  (8) 

which  is  the  probability  that  exactly  k  ones  occur  of  the  possible  m  ones 
present.  For  this  law,  the  mean  Is  mp  and  the  variance  is  mp(l-p).  Mow  each 
element  of  the  x  vector  is  formed  from  target  and  background  detections. 
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number  of  possible  false  alarm  or 
background  detections 


(10) 


We  will  use  the  notation  E(.)  as  the  expectation  operator.  The  variance 
Var(x)  is  defined  as 

Var{x)  =  E(x2)  -  E2(x)  vll) 

For  the  binomial  law 

E(V  =  tjPt  +  {n~V  pb  (12) 


and 

VaHxj)  .  tjPt  (1-Pt)  *  (n-tj)  ?b  (1-Pb)  (13) 
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so  that 

E(S)  .  S0Pt  *  (n2-S0)  Pb  (15) 

Var(3)  =,  S0Pt  (1-Pt)  *  (nZ-S0)  Pb  (1-Pb)  (16) 

(NOTE:  Although  there  is  partial  correlation  in  forming  the  sum,  independence 
i s  assumed. ) 

The  centroid  error,  js  defined  by 
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The  computation  of  E(x)  involves  the  ratio  of  two  random  variables  U  and  V, 
where 
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We  invoke  the  approximation, 
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From  Equations  (23) 
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By  independence, 
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This  result  indicates  that  when  a  centroid  error  results  (i.e.,  xq  is  not 
equal  to  zero),  the  process  produces  a  biased  answer  which  is  proportional  to 
the  error  and  the  probability  of  backgrouncTaetecti on . 

From  Equation  ( 34) , 


U0  =■  S0  <Pt-Pb>  x0 

Now,  to  simplify  examination  of  the  variances,  consider  the  50  percent 
detection  case  with  Gaussian  additive  noise. 
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Let  the  target  have  an  amplitude  A  above  the  background  and  the 
background  be  Nonna!  with  mean  4  and  variance  a2.  The  threshold  is  defined 
as 

T  =  u  +  A/2  (43) 

From  these  definitions. 
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From  Equation  (30) 
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Let  the  target  to  gate  area  be  given  by 
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For  the  tracking  case,  xq  driven  to  zero  so  that  the  residual  tracking 
noise  approaches, 

Ml-PJ 

Var(x)  i  - l - ± - 7  (64) 

12[l-Pt+R(2Pt-l)r 

Assuming  that  the  detection  Drobabi'.  ity  >  0.5,  the  form  of  the  equation 
indicates  that  the  variance  decreases  for  eitRer  an  increase  in  the  detection 
probability  or  an  increase  in  R.  Since  the  detection  probability  is  limited 
by  s-'gnal -to-noise  ratio,  the  result  suggests  that  the  gate  be  adjusted  to 
accommodat-*  the  target  and  that  the  best  result  is  obtained  when  R  is 
maximized. 

The  modelling  equations  are  summarized  below: 
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A  Monte  Carlo  experiment  was  conducted  using  a  200-sample  pattern  of 
8  by  3  squares  in  a  16  by  16  sample  matrix.  The  signal -to-noise  was 
established  at  the  values  shown  in  Table  I  and  the  statistics  were  compiled 
and  tabulated  (the  pattern  was  shifted  one  pixel  to  the  right  and  down). 


Using  the  measured  values  of  Pt  and  Pb  and  xq  =  1,  the  ^ 

Equations  (66)  through  (71)  were  used  to  compute  the  bias  error  E(e  x)  and 
the  standard  deviation  of  the  centroid  error. 
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These  data  are  shown  in  Table  II  and  indicate  excellent  agreement. 
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NPTS  = 2  0  0 
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3  2  A  1  ) 


:  '  .  F8  3  . 
1  .  F  8  3  , 

, F8  .  3  ) 

F  8  3  ,  4  X 


2  X 
2  X 


S I CMAX  = 
S I GMAY = 1 


'  PB=  '  ,  F8  3  ) 


SUBROUTINE  GNO I S E ( DC  .  S I GMA ,  I  NO  I SE ) 
DIMENSION  S  <  2  5  6  )  ,  I  NO  I S  E (  2  0  4  8  ) 
COMMON  PI 
P 1  =  1  / SORT  <  2  .  *  P I  > 

X  =  -  5 

H»  5 /SIGMA 
DO  10  I  =  1  ,  2  5  6 
Al=( X-DC ) / S I GMA 
A  2  =  (  X  f  5 -DC  >  /  5  I  GMA 
FlaEXPX— <Al*Al  )/2  > 

F2*EXP<— (AZ*AZ>  /  2  ) 

A  3  a (  X  + 1  -DC)  /SI GMA 
F3=EXP(~(A3*A3)/Z  ) 

3  «  I  /=H*Pl*(Fl+4  .  *F2tF3  )  /  3 
:<=Xt! 

Ml  =  1 

DO  20  1=1,256 

N  =  2  0  4  8  *  S  <  I  ) 

M2  =M  1  *•  N 

!  F  <  M2  GT  2  0  4  8  )  M2  =  2  0  4  8 


■  *•  -  ^ 


'  ,  F8  3  ) 
,f!  3) 


Ml 


3  0 

2  0 


DO  3  0  J  M  1  ,  M2 
I  NO  ISE(  J  '  I  -  1 
N  J  -M2  4-  1 
CUNT  I NUE 
R  E  TU  R  N 
END 

FUNCTION  I  GAUSS  '  I  NO ISE) 
DIMENSION  I  NO  I  SE  <  2  0  *1  6  ) 
M-  2  0  4  B  *  RAN  <1  )  +•  l 

ICAUSSalNOISEiM) 

RETURN 

END 
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13  D  I  3  P  T{  F  O  R 
P R O  G RAM  D I S  P R 

DIMENSION  IPAT;  H  i  ,  V,\LS  (  2  ,  2  ) 

D  I  MENS  I  ON  V  I  (  1  0  )  VMlQ) 

COMMON  PI 

DATA  VI/.  7’  45.  .823,  .887,  .922.  9  5  4,  971,  .987, 
DATA  V  2  >  18  4,  1  2  7,0  8  3,  0  4  2,0  2  1,  0  0  2  0  ,0 
DATA  IP AT /5*0, 8*1. 3*0/ 

P  I  =  4  .  *ATAN<  1  , 

SNR  = 1  ,  5 

DO  9  0  I SN= I  ,  9 

3 1  =  SNR /  2  . 

S  2  a  -  3  1 
PT=V 1  <  I  ON ) 

P  B  =  V  2  (  I SN ) 

DO  9  9  NCASE  =  1  ,  2 
M  -8 

I  F <  NC ASE  .  EQ .  2  )  H  =  0 
AM  =  1  6 
S=  6  4  . 

AMX=H*PT+<AM~H)  *PB 

SIG2X  =  H*PT*  (  1  ,  -PT)+<AM-H)*PB*(  1  -PB  > 

AMY  =  P  T  *  S  + ( AM  *  *  2  - S )  *PB 

3IG2Y  =  S*PTM  1  .  -PT)  +•  <  AM*  *2-S)  *PB* (  1  .  -PB) 

3  2  X  =  SORT ( 3  I G2  X ) 

S  2  Y  =  SORT  < S I G2 Y ) 

R  =  AMX / AMY 

R5S=R*SQRT<  3  I  G  2  X  /  (  AMX  *  AMX  )  t-3  I  G2  Y  /  (  AMY*  AMY  )  ) 

VALS (  1  , NC A  S  E )  =  R 

VA  L  S (  2  , NC  ASE)  =RSS*RSS 

CONTINUE 

SI  =0 

32  =  0  , 

DO  20  1=1,16 

A  =  F  LOAT  <  I  ) -8  .  5 

J=  1 

IF  <  I P  AT  <  I  )  .  EQ .  0  >  J  =  2 
3  1=31 +VALS (  1  , J  )  *  A 
S2=S2  +-VALS  (  2  ,  J  )  *  A  *  A 
3  1=31-1  . 

3  2  = SORT (32) 

WR ITE  (  2  ,  I  0  0  )  SNR.Sl  ,S2 
FORMAT (  I  X  .  3  F 6  3  ) 

SNR  a SNR  ,  5 

FORMAT (  1 H  t  ) 

WR I TE  <  2  ,  2  5  0  ) 

END 

FUNCTION  GAUSS ( X ) 

DOUBLE  PRECISION  B(5> ,P,T,T1 ,S 
COMMON  PI 

DATA  B/ .31938153, -.356563782,1  781477937 

,  -  1  321255978,1  33027  4-129/ 

P -  2316419 

Z=EXP<-(X*X>/2  )  / SORT (  2  * P  I  ) 

T -  1  /(I  +P  *DBLE ( ABS ( X  )  )  ) 


9  9  3 
.  0  / 


781477937 


T  1  -  T  1  *T 

0  SNC  L  <  S  )  *  E 

I  F  (  X  LT  0.)  G  =  1  -G 

GAUSS-  1  - G 

RETURN 

END 

SUBROUTINE  P  R  TV  A  L  <  L  V  A  L ) 

LOGICAL  LVAL 

I F ( ( INP ( X ' EE ’ ) . AND . 1 )  NE . 0 )  GOTO  10 
CALL  OUT ( X  '  EC  1  ,  LVAL ) 

CALL  OUT  C X  1  EF  '  ,14) 

CALL  OUT ( X  '  EF  1  ,  1  5  ) 

RETURN 
END 


A  .'Tv  PE  B  DISF’fi  FOR 

PROGRAM  DISPR 

D I  MENS  I  ON  I  PAT <  1 6  )  .  V A L S  <  2  ,  2  ) 

COMMON  PI 

DATA  [PAT / 5  *  0  .  8  *  1  ,  3  *  0  / 

P  I  =  9  *ATAN  (  1  .  ) 

SNR  =  1  .  5 

DO  98  I SN- 1 , 9 

3  1  -  S  N  P>.  /  2  . 

S  2  *  -  S  1 

PT  -  1  .  -GAUSS ( S2  ) 

P  B  =  1  .  -GAUSS ( S  l  > 

DO  9  9  NC  A  S  E  =  1  ,  2 

H  a  8  . 

I F ( NC AS  E  EQ .  2  )  H  =  0  . 

AM  =  1  6 
S=  6  4  . 

AMX=H*PT+( AM-H) *PB 

SIG2X=H*PT*  (  1  .  -PT)  +  UH-H  '  *PB*  (  1  -PE) 

AMY  =  FT* S  h  <  AM*  *  2 -S  >  *PB 

SIC2Y=S*FT*<1 ,-PT)+(AM**2-S)*Pa*U  -PB) 

S  2  X  =  SORT ( 5  I G2  X  ) 

S2Y  =>SQRT  (  SIG2Y) 

RrAMX / AMY 

RSS  =  R  *  SORT  <  SIG2X  /  (  AMX*AMX  )  +-SIG2Y  /  (  AMY*AMY>  ) 
VALS  (  1  ,  NCASE  )  --R 
V  A  L  S (  2  , NCASE)  -RSS* RSS 
”9  CONTINUE 

S  1  =0  . 

S2  =  0  . 

DO  20  1=1,16 

A= FLOAT! I ) -8  5 

J=  1 

I F  <  I P AT (  I  )  EQ .  0  )  J  =  2 
S1=S1+VAL3< 1 ,J)*A 
20  S2=S2+VALS<2,J>*A*A 

S  1  =S  1  -  1 

S  2  =  SORT ( S2-S1 *S1  ) 

WR ITE  <  2  ,  10  0)  SNR  .SI  ,  S  2 
1  0  0  FORMAT!  1 X  ,  3F6  3  ) 

98  SNR = SNR f  5 

?T0  FORMAT  (  1  H  1  ) 

WR I TE (  2  ,  2  5  0  ) 

END 

r 

FUNCT ION  GAUSS  <  X  ) 

DOUBLE  PRECISION  E(5) ,P.T,T1 .3 
COMMON  PI 

DATA  B/.  31  >  38153,-  3  56563782,1  781-3  77937 

1  ,  -  1  821255978,133027  3  329/ 

P  =  2  3  1  6  9  1  9 

Z=EXP(-(X*X)/2  ) /SORT! 2  *PM 

T  =  1  /  (  !  +P*DBLE(ABS<X)  >  ) 

T  1  =  T 
S  =  o 

DO  : 0  1=1,5 

S  =  S  T  1  *  B  <  I  . 

10  Tl-Tl  *  T 

G  =  SNG  L  '  3  )  *  7 


IF  <  X  LT  0  )  Q-=t  -O 

' A  U  S  5=1  -  a 

■IE  TURN 
END 


SUBROUTINE  P MTV A  L ( LVAL  ) 

LOGICAL  LVAL 

I  F  (  (  I  NP  <  X  ’  EE  ’  )  .  AND  1  !■  ME  0  >  GOTO 
CALL  O  UT (  X  1  E  C  1  , LVAL  ) 

CALL  OUT ( X  1  EF  1  .10) 

CALL  OUT ( X  1  EF  1  ,15) 

RETURN 

END 


^  i-v  -w'.  ..  '  ± *1— .  *w»:. 


BINARY  CORRELATION  TRACKER 


The  correlation  tracker  to  be  analyzed  below  assumes  the  following  form. 

The  tracker  is  a  "searching  type"  wherein  a  reference  is  taken  in  one 
video  frame  and  a  picture  sample  is  search  in  subsequent  frames  to  determine  a 
location  where  some  sub-matrix  in  the  picture  sample  best  matches  or  maximally 
correlates  to  the  reference.  Errors  are  determined  by  noting’ where  in  the 
raster  the  picture  sample  is  taken  and  where  in  the  sample  the  best  match 
occurs.  The  above  description  is  served  by  an  example  which  is  used  in 
subsequent  experiments. 

Suppose  the  picture  sample  is  a  16  by  16  binary  matrix  and  tnat  at  the 
initializing  time,  the  reference  is  a  centrally  located  12  by  12  matrix  within 
the  picture  sample.  The  matrices  are  shown  -*0  Figure  1.  In  subsequent  frames 


s 


the  reference  matrix  is  placed  in  the  upper  left  hand  corner  of  the  picture 
sample,  performs  a  correlation  and  is  then  moved  one  pixel  to  the  right  for 
the  next  correlation,  etc.  as  shown  in  Figure  2. 


In  all,  25  positions  of  the  reference  will  cover  the  picture  sample.  The 
resultant  correlations  are  given  in  a  5  by  5  matrix  C.  The  boresight  position 

corresponds  to  C33.  Suppose  maximum  correlation  occurs  at  Cij.  The 
errors  are  then  given  by 


j  - 


o 

sj 


A 


3  -  i 


(1) 


Figure  2. 


(4) 


2  r  2 

J  hi to vi 


Now  suppose  that  m  of  the  elements  of  u  match  m  of  the  elements  of  v  one  for 
one  with  m  £  n.  Then, 

n 

T  u.v,  =  ,n  -  (n-m)  ( 

1=1  1  f 

That  is,  the  remaining  n-m  elements  have  opposite  sign.  Consequently, 


)  ( u ,  v )  = 


2m  -  n 


If  all  of  the  elements  match,  the  correlation  is  one  and  if  none  of  the 
elements  match,  the  correlation  is  minus  one.  Now,  working  with  conventional 
binary  numbers;  i.e., 

V  v1  €  j  0,1  } 

t  Vi  s  n  “  2  Z  fi  ©  vi  (i 

i-1  1  1  i -1  1  1 

where  the  symbol  (£)  implies  the  "exclusive  or"  function.  Then, 

2  *  number  of  "matches"  ,  /, 

p(u,v)  =. - - .  I 


The  "exclusive  or"  truth  table  is  as  follows: 


U  ©  V 


Now  given  the  matrix  R,  we  form  a  reference  vector  r  as  follows.  Let 


k  =  j  +  ( i  - 1 )  n 
and 

=  C-j  ]  i  ,j=l,  • « •  ,n 
;  k=l , 2 , . . . , n  2 

Similarly,  define  a  sub-matrix  of  S  as  follows 

yk^  _  c 

ij  "  t+k-l,jt*-l 

i  ,j=l,2, . . .  ,16 
k ,  j? = 1 , 2 , .  ,  5 


Then 

t  (I-  /]  _  rk-^  c 

V  '  'ij  "  J1+k-ltj+i-l 

where 

m  =  j  +  12(i-l) 

Then 

ck£  =  p  (  r,t(k,£)  j 


(10) 

(11) 

(12) 

(13) 

(14) 

(15) 


The  best  correlation  implies  that 

ck*£*>ck^  for  all  kj=l,2s...,5  (16) 

and  (k*,^*)  is  the  track  point. 

The  basic  correlation  mechanism  is  rather  straightforward;  the  analysis, 
however,  is  quite  tedious.  The  search  process  produces  a  decision  tree  which 
must  be  evaluated  to  determine  the  estimated  best  point.  As  will  be  shown, 
the  best  point  determination  is  biased  according  to  the  order  in  which  the 
points  are  searched. 

For  the  example  given,  25  points  must  be  searched.  Let 

Pk  =  Prob  {position  k  is  chosen  |  ;  k=l,...,25  (11) 


<34 


Further,  let 

lUm.n)  =  Prob  j  correlate  on  at  position  m 

>  correlation  at  position  n|  (18) 

At  the  start,  of  the  process, 

P,  =  R( 1,2) 

1  (19) 

P2  =  R{ 2, 1) 
and 

R( 2, 1)  =  1  -  R ( 1,2)  (20) 

with  n  =  2 

The  conditional  probabilities  can  be  computed  b.y  the  following  iteration: 


Vi  ■  i  R(nn>i,pi 

Pk  =  (Ukmn+DP^  ;  k=l,...,n 

n  =  n+l 

Equation  (21)  is  interpreted  as  follows: 


(21) 

(22) 

(23) 


The  n+l  state  is  chosen  provided  that  the  probability  of  the  n+l  state  is 
greater  than  some  state  i  given  that  the  state  1  had  been  previously  chosen. 
Equation  (22)  implies  that  state  k  is  chosen  provided  it  indicates  a  greater 
probability  than  the  current  state  and  the  probability  that  it  had  been 
previously  chosen  as  the  optimum  state.  The  process  stops  when  n+l=m  (or  25 
in  the  current  example).  At  any  stage  j  in  the  process,  the  probability  of 

choosing  one  of  the  states,  Py  is  given  by 


pJ 

T 


=  y 

P. 

i 

■L 

R(j 

1 

_  \ 

P.  L 

i4l 

i 

Pi  =  t  +  t  RCi  .J-DP,-  ; 

T  i=l  1  i=l  1 


;  j  >  2 


:(i,j-i)J  =  f  p, 

iTn  1 


(24) 


(25) 


(26) 


Prob  { 1  chosen  [  =  Pj. 


R14R13R12 


Prob  |  2  chosen  } 
Prob  {  3  chosen  | 
Prob  | 4  chosen  \ 


P2  3  R21R23R24 

P3  3  R34(R31R12*R32R21) 

P4  s  R41R13R12  +  R43(R31R12+R32R21)  +  R42R23R21 


PT  .  Pl>P2  P3>P4 

PT  -  R12(«14R13+R34R31+R41R13+R43R31) 

+  R2l^R24R23+R34R32+R43R32+R42R23^ 


PT  =  RI2  R13^14+R41^+R31^R34^R43^ 

*  R?1  R23^R24fR42^+R32^R34+R43^ 

(34) 

PT  =  R12(R1;^31)  +  R21fR23+R32) 

(35) 

P.j.  a  R12  +  R21  "  ^ 

(36) 

that  to  search  out  this  tree, 

aLsAl  ,  325 

probabilities  have  to  be  determined  (although  some  symmetries  reduce  this 
number) . 


The  problem  is  to  determine  the  probability  that  one  position  is  better 
than  another;  however,  this  problem  is  cor.vpl (sated  oy  the  fact  that  the  data 
in  nearby  positions  is  highly  correlated.  The  algorithm  >jsed  to  compute  the 
probabilities  is  shown  in  FORTRAN  code  in  Figure  4,  i.e„,  SUBROUTINE  E'/AL. 
The  parameters  kl,  k2,  k3,  and  k4  are  reference  position  offsets.  A  diagram 
of  the  matrix  positions  is  shown  in  Figure  5. 


The  algorithm  computes  the  difference  between  the  correlation  values 
produced  by  the  two  reference  positions  (AM)  and  the  variance  of  the 
difference  measure  (VAR).  The  probability  that  one  position  is  better  than 
another  is  the  probability  that  the  normalized  difference  value  is  greater 
than  zero.  The  assumption  is  that 


P(x) 


1 


(x-AMr 

~mr~ 


/2ir  VAR 


,-OG 


PROB  {x  >  0} 


J 

o 


(x-AM)2 
2  VAR 


dx 


oo 


u2 

~1 


-  AM 

sm 


(37) 


(38) 


PROB  {  x  >  0  } 


1 

/27 


e 


du 


(39) 


IPII.1) 

"TaEUEIS 


1 R  E  F  ( 1 2 .  J  2 ) 

1 

Figure  5. 


and  for 


(40) 


(41) 


The  results  of  these  computations  are  shown  in  Table  III.  The  striking  result 
is  that  the  errors  are  biased  according  to  the  path  taken  by  the  search 
pattern;  i.e.,  to  the  right  and  down. 

To  test  the  theory,  a  Monte  Ca^o  test  was  run  using  200  samel es  at 
various  SNR.  The  results  are  shown  in  the  following  figures.  The  binary 
pattern  shown  corresponds  to  the  initial  reference.  Data  indicate  that  the 
errors  are  indeed  biased.  A  5  by  5  error  matrix  is  also  shown  which  indicates 
the  number  of  times  over  the  200-sample  run  that  a  particular  error  position 
occurred  (when  divided  by  200,  this  matrix  approximates  the  probability 
matrices  shown  in  Table  111).  Although  the  analysis  is  somewhat  pessimistic, 
tii a  agreement  is  good  and  the  implied  data  trends  are  apparent..  Another 
interesting  result  is  that  for  rather  moderate  SNR  (  *  2.5),  the  correlator 
essentially  makes  no  error  and  can  always  correctly  identify  the  correct 
target  position;  a  significant  error  is  made  with  the  centroid  measure  for  the 
same  SNR. 


<3?(x) 

PROB  |  x  >  0  | 


/75r  '•'-GO 


u2 
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e  du 
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:  wrfjs  :  ou  i>c.  y.oa  e  ■  :■  5  •  ..r  l  ;  =0  ■  2 5  ■ 

D  L  MEM  £  1  JM  OR  OB  •  2  '.■  .  i.  .  '  f»OS  1  2  2  .  2 

;3  ;;|rvj  r  q.j  r  i>  <  i  *.  t  .  •. 

:  .  -tMOM  ?5'i 

MATA  I  IT  t  0  ,  :}  0  ,  3  .  3  .2.3  0  3.0.  5  S  3  .0.0  -  - 

i  .3.0.5.  !  0  t  *.  2  0  .  0  / 

3  '  |  ‘  AT  A  M  '  1 

o  a  o  v  1  '  £  0  p  2  cs  !  i 

D  O  5  0  0  i  i  ,  16 

0  0  3  0  0  J=»  ;  ,  1  6 

■  F  C  I  .  J  >  a  C 

•  "  ■'  i  I  C  S  .  4  )  AMD  .  i  I  LE  13'  AMD  (JOE  j  >  AMD  >  J  L  E 

r  p  ■:  j  .  j  '  -  i 

GO MT I NUE 
DO  510  1=3. Id 

0  0  5  1  0  J  =  3.1-1 

i  r  e  r  <  i  - :  ,  j  -  2  >  =■  r  p  i  ,  ■  7  > 

FORMAT i / 1  > 

FORMAT  1  H  1  .' 

3HR= ! . 0 

DO  10  I  PAR  -  1  ,5 
S  I  GMA  =  6  4  .  /  SNR 
DC =3  ‘SIGMA 
S  !  GMA  L  =  <  SNR  t-  3  >  '  3  GMA 

I  T s  (  DC  ~S  I  GN'AL  >  /  2 

CALL  GNO I SE -  DC  . S  t  GMA . FT . PE  .  IT)' 

VR ITE ( 2  ,  1 00  >  SNR .  PT .  P9 


DO  5  7  5 

i  =  i 

,  2  5 

PRO  3  4  I 

i  >  = 

0  . 

DO  5  7  3 

i  - 1 

c. 

DO  5  7  6 

j--  i 

c- 

-  J 

M=  J  *  5  *  i. 

i  - 1 

\ 

IPOS  M'N 

i  ’  - 

l 

IPOS ' N , 

2  5  •= 

j 

DO  5  8  Q 

I  ! 

.  2  <\ 

I  1  —  I  j-  L 

do  sac 

J  -  I 

1  .  2 

K  1  -  I  P  0  S  -  I  .  1  > 

K  L  -  I  P  0  £  1  [  ,  2  ) 

J J  =  J L I  ST < J  ' 

I  F  -  (  J  J  1-0  0)  OR  (JO  LT  I!  GOTO  573 
?  0  -  P  R  O  3  (  ■  .1  J  > 
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